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Abstract. We define Szcgd coordinates on a finitely connected smoothly 
bounded planar domain which effect a holomorphic change of coordinates on 
the domain that can be as close to the identity as desired and which convert 
the domain to a quadrature domain with respect to boundary arc length. 
When these Szego coordinates coincide with Bergman coordinates, the result 
is a double quadrature domain with respect to both area and arc length. 
We enumerate a host of interesting and useful properties that such double 
quadrature domains possess, and we show that such domains are in fact dense 
in the realm of bounded C°° -smooth finitely connected domains. 



1. Introduction 

The unit disc -Di(O) is the quintessential example of a double quadrature do- 
main. Holomorphic functions h on the disc satisfy two quadrature identities: 

n h{0) =11 hdA and 2nh{0) = / h{z) ds, 

JJd 1 (0) A,Di(0) 

whenever these integrals make sense. We will show in this paper that double 
quadrature domains are rather commonplace among smoothly bounded domains 
in the plane. 

A bounded domain Q in the complex plane is a quadrature domain with respect 
to area measure if there exist finitely many points {wj}j =1 in the domain and 
non-negative integers rij such that complex numbers Cjk exist satisfying 

j=l k=0 

for every function h in the Bergman space H 2 (Q) of square integrable holomor- 
phic functions on Q. Here, dA denotes Lebesgue area measure. We shall say that 
such a domain is an area quadrature domain. Aharonov and Shapiro [1] proved 
that such domains have piecewise real-analytic boundaries (see also |14j). 

A bounded domain Q bounded by finitely many non-intersecting piecewise 
C 1 -smooth Jordan curves is a quadrature domain with respect to boundary arc 
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length measure if there exist finitely many points {wj}^ =1 in the domain and 
non-negative integers rij such that complex numbers Cjk exist satisfying 

„ N rij 

/ hds = ^2^2c jk h^{ Wj ) 

for every function h in the Hardy space H 2 (bfl) of holomorphic functions with L 2 
boundary values. Here, ds denotes the boundary arc length measure. We shall 
say that such a domain is a boundary arc length quadrature domain. Boundary 
arc length quadrature domains have C^-smooth real-analytic boundaries (see 
[15], where it is also shown that the boundary regularity assumptions can be 
greatly relaxed). 

After some pioneering work of Avci [2] and Shapiro and Ullemar [21], the sec- 
ond author [15] studied planar quadrature domains with respect to boundary 
arc length from the point of view of half-order differentials and Riemann sur- 
face theory; he characterized such domains and showed that they are dense in 
the category of bounded domains in the plane bounded by finitely many non- 
intersecting Jordan curves. Here, we shall reframe these results in more elemen- 
tary terms and explain how to view boundary arc length quadrature domains in 
terms of Szego coordinates, which are the analogue of Bergman coordinates as 
developed in [HI QUI EH] • Because our results are most interesting in the category 
of C°°-smoothly bounded domains, we shall take the royal road and restrict our 
attention to domains of this form. 

For more information about the state of research on quadrature domains, see 
[T7] and the volume JT3] that holds it, and [20] . We also remark that it is possible 
to study quadrature domains with respect to harmonic functions in the plane, 
and to prove density results in this context (see [T5] and Sakai [T9]). 

Suppose that fl is a bounded domain bounded by n non-intersecting C°°- 
smooth Jordan curves. Let A^iVt) denote the subspace of C°°(Q) consisting of 
holomorphic functions. Let K(z, w) denote the Bergman kernel associated to Q 
and let S(z,w) denote the Szego kernel. Let K°(z,w) also denote K(z,w) and 
let K m (z,w) = -^^K{z,w). As in [IT] , we define the Bergman span associated 
to Q to be the complex linear span /C of all functions h(z) of z of the form 
K m (z, a) as a ranges over Q and m ranges over all non- negative integers. If U 
is an open subset of Q, let K,u denote the complex linear span of functions of z 
from {K(z,a) : a e U}, and given a point a in Q, let K. a denote the complex 
linear span of functions of z from {K m (z, a) : m = 0, 1, 2, . . . }. Similarly, define 
the Szego span to be the complex linear span S of all functions h(z) of z of the 
form S m (z,a) as a ranges over Q and m ranges over all non-negative integers, 
and define Su and S a in the analogous way. 

We now list two theorems that will serve as the heart of the paper. Note that 
if / : Q\ — > Q2 is a biholomorphic mapping between finitely connected domains 
in the plane, then it is well known that there is a holomorphic branch of ^/J , on 
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Theorem 1.1. Suppose {l\ and f2 2 are bounded domains in the plane bounded by 
finitely many non-intersecting C 00 -smooth Jordan curves. Suppose further that 
f : Qi — > Q2 is a biholomorphic mapping. Then Q2 is a boundary arc length 
quadrature domain if and only if f is equal to the square of an element of the 
Szego span associated to fit, or equivalently, if and only if a/7 7 is in the Szego 
span. 

We shall call a biholomorphic mapping / as given in Theorem 11.11 such that 
a/J 7 is in the Szego span a Szego coordinate. This is a direct analogue of the 
definition of Bergman coordinate in [TT] where a biholomorphic map / is called 
a Bergman coordinate if and only if /' is in the Bergman span. 

The next theorem allows us to use Theorem 11.11 to approximate domains by 
boundary arc length quadrature domains. 

Theorem 1.2. Suppose that Q is a bounded domain bounded by finitely many 
non-intersecting C°°-smooth Jordan curves. Then the Szego span S associated to 
Q is dense in A°°(Q). Furthermore, Su is dense in A°°(Q) for any open subset 
U ofQ, and S a is dense in A°°(Q) for any point a £ Q. 

Theorem 11.21 was proved in [I] . 

Theorems 1 1 . 1 1 and 1 1 . 2 1 and results of [TT] and [TH] can be combined to obtain the 
following results. The next theorem shows that boundary arc length quadrature 
domains are dense in the realm of bounded finitely connected smoothly bounded 
planar domains. 

Theorem 1.3. Suppose Q is a bounded domain bounded by finitely many non- 
intersecting C°° -smooth Jordan curves. There is a biholomorphic mapping f 
which is a Szego coordinate as close to the identity in C°°(fi) as we please. 
Hence, /(O) is a boundary arc length quadrature domain which is C°° close to 
Q. 

The proof of Theorem 11.31 was sketched in [18, p. 285-6]. We shall flesh out the 
argument here in §H to make this paper more comprehensible and because we 
will need the fine details in order to generalize the result to double quadrature 
domains. 

Theorem 1.4. Suppose Q is a bounded n-connected domain bounded by n non- 
intersecting C°° -smooth Jordan curves. There is a biholomorphic mapping f 
which is both a Szego coordinate and a Bergman coordinate as close to the identity 
in C°°(fl) as we please. Hence, /(fi) is both a boundary arc length quadrature 
domain and an area quadrature domain, i.e., a double quadrature domain, which 
is C°° close to n. 

We will prove Theorem 11.41 in §3] in case the domain is simply connected and 
prove the general case in §SJ 

The Riemann mapping theorem says that every simply connected domain can 
be mapped to a special double quadrature domain, the unit disc. Our results can 
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be thought of as a generalization, or maybe even an improvement. Every finitely 
connected domain in the plane can be mapped to a nearby double quadrature 
domain. 

Double quadrature domains have many remarkable properties in common with 
the unit disc, some of which we enumerate here. 

• The boundary is given by finitely many non-intersecting real algebraic curves 
which are C°°-smooth and real analytic. (It is shown in [14J that the boundary 
is essentially given by the zero set of an irreducible polynomial P(z,w) on C 2 . 
In fact, the boundary is equal to {z : P(z, z) = 0} minus perhaps finitely many 
points, and P(z,w) must have a a certain special form.) 

• The complex polynomials belong to the Bergman span. 

• The complex polynomials belong to the Szego span. (In fact, this property 
characterizes double quadrature domains. See Theorem 16. 21 ) 

• The complex unit tangent vector function T(z) is a rational function of z and 
z for z in the boundary (which of course is equivalent to being a rational function 
of Re z and Im z). The function T(z) is also the restriction to the boundary of a 
meromorphic function on the double without poles on the boundary. This means 
that T(z) can be extended to the domain as either a meromorphic function with 
no poles on the boundary, or an anti-meromorphic function with no poles on 
the boundary. (A more precise statement of this known fact can be found in 
Theorem 16.11 ) 

• The Schwarz function S(z), which exists and extends to be analytic on a 
neighborhood of the boundary by virtue of the real analyticity of the boundary, 
also extends to be meromorphic on the domain. It has no poles on the boundary, 
and it is an algebraic function. Because S(z) = z on the boundary, both z and 
S(z) extend meromorphically to the double of the domain. It was noted in [T3] 
that all the meromorphic functions on the double are generated by the extensions 
of z and S(z) to the double. 

• Both the Bergman kernel K(z, w) and the Szego kernel S(z, w) are rational 
functions of z, S(z), w, and S(w). Hence, they both extend to the double cross 
the double as functions which are meromorphic in z and anti-meromorphic in 
w. If the domain is multiply connected, the two kernels are rational functions of 
fi( z ), f2( z ), fi(w), and f2(w), where fi and f2 are two Ahlfors maps associated 
to two generic points in the domain. In the simply connected case, the kernels 
are even more elementary (see below). 

• Both the Bergman and Szego kernels are rational functions of z, z, w, and 
w when z and w are restricted to the boundary. 

• The Kerzman-Stein kernel A(z, w) is a rational function of z, z, w, and w 
for z and to on the boundary. 

If the domain is a simply connected double quadrature domain, then even 
more can be said. Let / denote a Riemann map associated to a point in the 
domain (which maps the domain one-to-one onto the unit disc). 
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• / extends to the double because f(z) = 1/ f(z) on the boundary. Hence, 
it follows that f(z) is a rational function of z and S(z), and so it is a rational 
function of z and z on the boundary. It is also algebraic. Since the extension 
of / to the double generates the meromorphic functions on the double, it also 
follows that z and S(z) are rational functions of f(z). (This is another way to 
deduce that / has a rational inverse, a well know fact about simply connected 
quadrature domains with respect to 

• The Bergman kernel K(z, w) and the Szego kernel S(z, w) are rational func- 
tions of f(z) and f(w). 

• The Poisson kernel p(z, w) is the real part of a function that is rational in z 
and S(z), and w and w (for z in the domain and w on the boundary). It is also 
the real part of a rational function of f(z) and w and w. It is even the real part 
of a rational function of f(z) and f(w). 

We shall explain these results in fj6l In £0 we reveal a connection between 
the Schwarz function, quadrature domains, and proper holomorphic mappings 
onto the unit disc. We prove, for example, that a domain is an area quadrature 
domain if and only if, on the boundary, z/z is equal to the boundary values of 
a quotient of proper holomorphic mappings to the unit disc (Theorem 17.11) . A 
smooth domain is an arc length quadrature domain if and only if the complex 
unit tangent vector function on the boundary is equal to the boundary values of 
a quotient of proper holomorphic mappings to the unit disc (Theorem 17. 2p . 



Suppose that / : Vt\ — > £l 2 is a biholomorphic mapping between bounded 
domains in the plane bounded by finitely many non-intersecting C°°-smooth 
Jordan curves. Let F = f" 1 . It is well known that a holomorphic square root 
of /' exists on Qi and a square root of F' exists on Q 2 - Since /' and F' extend 
C°°-smoothly to the respective boundaries and are non- vanishing there, the same 
is true of their square roots. Notice that the change of variables formula yields 
that 



denote the Hardy space inner product on Qj, j = 1,2. The last identity can be 
written 



2. Proof of Theorem 11.11 




when u is in the Hardy space H 2 



(6ft 2 ). Let 
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To prove the first half of Theorem 11.14 assume that a/7 7 is m the Szego span 
associated to f^. If h is an element of the Hardy space H 2 (6O2), then 

/ hds= [ \f'\(hof)d8=(y/r{hof),y/f% 

and if \/J 1 is in the Szego span, then this last inner product yields a finite linear 
combination of values of \fj r {h o /) and its derivatives at finitely many points in 
Qi, which reduces to a finite linear combination of values of h and its derivatives 
at finitely many points in Q 2 . This shows that Q 2 is a boundary arc length 
quadrature domain. 

To prove the converse, suppose that fl 2 is a boundary arc length quadrature 
domain. Then, given any g G H 2 (bQi), 

(9, v r F)i = (VF^(goF),l) 2 , 

and since fl 2 is a boundary arc length quadrature domain, this last inner product 
yields a finite linear combination of values of \/F'(g o F) and its derivatives at 
finitely many points in Q 2 , which reduces to a finite linear combination of values 
of g and its derivatives at finitely many points in Q\. Thus, a/7 7 has the same 
effect as an element of the Szego span when paired against a function g in the 
Hardy space. It follows that y/J 7 must in fact be equal to that element in the 
Szego span. 

3. Density of quadrature domains in the simply connected case 

In this section, we give two proofs that boundary arc length quadrature do- 
mains are dense in the realm of bounded C°°-smoothly bounded simply con- 
nected domains, and that double quadrature domains are dense as well. The 
first proof introduces techniques that we will generalize to the multiply con- 
nected setting. The second proof is so easy, natural, and short that one could 
get the wrong impression that the multiply connected generalization should be 
easy and straightforward. 

Suppose that Q is a bounded simply connected domain in the plane bounded 
by a C°°-smooth Jordan curve. We may now construct a boundary arc length 
quadrature domain close to Q as follows. By Theorem 11.21 there exists a function 
h in the Szego span that is close to the constant function 1 in A°°(Q). (Of course, 
we take h close enough to 1 so as to be no n- vanishing.) Let if be a holomorphic 
square root of h on Q, and let / be a complex antiderivative of H, where we 
choose the constant of integration so that f(z) is close to the identity in A°°(Q). 
By choosing h close enough to 1 in A°°(Q), we may guarantee that / is close 
enough to the identity that / is one-to-one on Q and f(Q) is as C°°-close as we 
desire to Q. 

Now according to Theorem 11.11 /(fi) is a boundary arc length quadrature 
domain that is C°° close to Q. 
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We can repeat the above argument using S a for a point a in f2 instead of the 
complete Szego span. We obtain a biholomorphic mapping / to an arc length 
quadrature domain close to Q satisfying an identity, 

2 



A? 



/'= nr Cm s m (z,a) 



vm=0 

We want to see that /' is also in the Bergman span so that f(Q) is also an area 
quadrature domain via Theorem 1.3 of [IT] . 

Let L(z,w) denote the Garabedian kernel associated to Q and let L m (z,w) = 
(d/dw) m L(z,w). Let A(z,w) denote the Schiffer kernel function, which is mero- 
morphic in z and w with a double pole at z = w and which is related to the 
Bergman kernel K(z, w) of Q via 



K(w,z)T(z) = -A(w,z)T(z) 

when w £ Q and z G bfl. Here T(z) denotes the complex unit tangent vec- 
tor function pointing in the direction of the standard orientation. The ba- 
sic properties of these kernel functions are described in [5]. We also define 
A m (z,w) = (d/dw) m A(z,w, 
determined later. The relationship 



Let bj denote complex coefficients that will be 



(3-1) 



S(w,z) = -L(z,w)T(z) 



for w G Q and z G bQ together with the relationship between K(z, w) and A(z, w) 
allow us to state that H(z)T(z) = —G(z)T(z) for z in bQ, where 



and 



H(z) 



G{z) 




M 



l (z,a) 



N 
■.m=0 



[z, a 



It can be deduced that ( J2f=o c jL 3 (z, a 




T r 



[z, a 



dz 



m=0 




z, a) 



is residue free since 



Or, 



lz,a 



dz, 



and this last integral is zero by Cauchy's Theorem. Since A(z,a) has principal 
part (— 1/tt)(z — a)~ 2 at z = a, it follows that a positive integer M and constants 
bj can be chosen so that G(z) has a removable singularity at a. Now both 
H(z) and G(z) are holomorphic on Q and extend smoothly to the boundary in 
such a way that H(z)T(z) = —G(z)T(z) on the boundary. We now claim that 
this can only be true if both H and G are identically zero. Indeed, functions 
of the form G(z)T(z) are orthogonal to holomorphic functions in the Hardy 
space and functions of the form H(z)T(z) are orthogonal to anti-holomorphic 
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functions. Thus, H(z)T(z) is orthogonal to holomorphic functions and anti- 
holomorphic functions. But every C°°-smooth function on bfl is the restriction 
to the boundary of a harmonic function in C°°(Q), and every such harmonic 
function can be written as h + h where h G A°°(Q). Since smooth functions are 
dense in L 2 (6fi), it follows that H = and G = 0. We can now state that 



and this shows that /' is in the Bergman span. Hence / is both a Szego and a 
Bergman coordinate and /(fi) is a double quadrature domain. 

We remark here that, after finding the proof above, we discovered that Avci 
[2] showed in his unpublished Stanford PhD thesis using other methods that 
an arc length quadrature domain that satisfies a one point quadrature identity 
would also have to be an area quadrature domain. We include our proof here 
because it shows the close and explicit connection between Szego and Bergman 
coordinates, and because we will improve upon it in £j5]when we generalize our 
results to smooth multiply connected domains. 

We now turn to the simple proof of density of double quadrature domains in 
the simply connected setting promised at the beginning of this section. The proof 
is based on a theorem of Aharonov and Shapiro for area quadrature domains and 
a theorem of Shapiro and Ullemar for arc length quadrature domains. Let Q be 
a bounded simply connected domain bounded by a C°°-smooth Jordan curve 
and let / denote a Riemann mapping of Q to the unit disc. Let F = Note 
that a branch of \ ~F' can be chosen which is analytic on the unit disc and which 
extends to be in A°°(Di(0)). Aharonov and Shapiro pQ proved that Q is an area 
quadrature domain if and only if F is rational. Shapiro and Ullemar [2TJ proved 
that Q is an arc length quadrature domain if and only if F' is the square of a 
rational function. To construct a double quadrature domain that is C°° close to 
fi, we may approximate y/~F' in A°°(D±(0)) by a polynomial P(z). Let p(z) be (a 
polynomial) antiderivative of P{z) 2 where the constant of integration is chosen 
so that p(z) is close to F' in A°°(Di(0)). If the first polynomial approximation 
P(z) is sufficiently close to vF', then the image of the unit disc under p(z) will 
be a double quadrature domain that is C°° close to Q. 

4. Density of boundary arc length quadrature domains in the 



In this section, we show how the argument in [151 P- 77-78] using half-order 
differentials and a special Runge Theorem on Riemann surfaces can be modified 
along the lines of the previous section to prove C°° density of boundary arc length 
quadrature domains in the category of smooth finitely connected domains. 

Let Q be a bounded domain bounded by n non-intersecting C°°-smooth Jordan 
curves. There is a function h in the Szego span that is as close to one in C°°(Q) as 
we desire. We will need the rational functions Rj(z), j = 1, . . . , n — 1 constructed 
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in [T51 p. 78] that satisfy the following special orthogonality conditions. Let jj, 
j = 0, . . . , n — 1 denote the boundary curves of Q where 70 denotes the outer 
boundary Choose a point Zj inside jj, j = 1, ... ,n — 1 (i.e., choose one point 
from each of the holes in Q). The rational functions are given by 

Rs{z) = 2m(! - Zj ) + {Z ~ Zj)Qj(zl 
where the Qj are any polynomials satisfying 

1 



Qj(zk) 



2ni(z k - zj) 2 

for k 7^ j (making Rj(zk) = for k ^ j). The functions Rj have the property 
that 



Rj dz = 5k j and / RjR m dz = 
'ik J"ik 
for j, k, m = 1, . . . n — 1. Since each rational function Rj is in A°°(Q), we may 
approximate it in A°°(Q) by a function Tj in the Szego span. Notice that it easy 
to find coefficients Aj such that the periods of (h — Y^j=i ^i-^i) 2 vanish because 
the system is 



f h 2 dz + 2^Aj f hRj dz = 



'Ik j = \ J lk 

for k = 1, . . . ,n — 1. Since h is close to one, the coefficient matrix 



hRj dz 

'ik 

is close to the identity and f h 2 dz is close to zero. Thus, the Aj are uniquely 
determined and are close to zero. If we modify the problem of making the periods 
vanish by replacing the Rj by our functions Tj in the Szego span, then the system 
becomes 

» re— 1 » n—l „ 

/ h 2 dz + 2 Aj / hrj dz + ^ AjA m I Tjr m dz = 0, 

Jlk j= i Jlk j tTn= x Jlk 

and since the coefficients in the linear part are close to the identity matrix and the 
coefficients of the quadratic terms are small, the implicit function theorem yields 
that the Aj are uniquely determined and are close to the small solution to the 
system with Rj in place of rj. Notice that by taking h to be sufficiently close to 
one, the coefficients A,- can be made small enough that h - A 

jTj is as close 



to one as desired. Thus, this shows that there is a function F = h — Y^j-i A j r j 
in the Szego span which is close to one whose square F 2 has vanishing periods. 
Now let / be an antiderivative of F 2 . Since F 2 is close to one, we may choose 
the constant of integration so that / is close to the identity. If h is close enough 
to one, then / will be close enough to the identity to make / biholomorphic on 
fl and Theorem 11.11 yields that / maps Q to a boundary arc length quadrature 
domain that is C°° close to Q. 
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5. Density of double quadrature domains in the general case 

We continue with the same domain and setup as in the previous section. In 
particular, 7^, j — 0, . . . , n— 1 denote the boundary curves of fl where 70 denotes 
the outer boundary. Let p — n — 1, and let fl denote the double of fl. We are 
now thinking of the curves jj as being curves on fl. Define p more curves 7^, 
j = p + 1, . . . , 2p, on fl as follows. Let Tj be a simple smooth curve that starts 
at a point on jj, enters fl, and terminates at a point on the outer boundary 70. 
We may choose the Tj so that no two of them intersect. In this case, the domain 
U given by fl — U^ =1 Tj is a simply connected domain. Let Tj denote the curve 
which is the image of the reverse of Tj in the reflected copy of fl in fl, and let 
7 P+:) be the closed curve in fl obtained by following Tj in fl and then Tj in the 
reflected copy of fl back to the starting point. In this way, we obtain a standard 
homology basis, jj, j — 1, . . . , 2p, for fl where the first p curves are the standard 
homology basis for fl and the last p curves wrap around the "handles" of fl and 
complete the homology basis for fl to a homology basis for fl. 

We will need the concept of a half-order differential on the double as used in 
[13] . In the present context, we will describe a meromorphic half-order differential 
h on fl as a pair of meromorphic functions (hi, h 2 ) where the two functions are 
meromorphic on fl, extend continuously to the boundary (and are allowed to 
take the value 00 at a boundary point), and satisfy 

(5.1) h x (z)T{z) = h&j 

on the boundary of fl. The symbol h stands for the pair (hi, h 2 )- These objects 
have the virtue that identity ( 15. ip allows the meromorphic differential h\ 2 dz on 

fl to be extended to a meromorphic differential on fl (via extension by h 2 dz on 
the anti-holomorphic reflected side). We let h 2 dz denote the extension of h\dz 
to the double. Another important feature is that if h = (hi, h 2 ) and g = (gi,g 2 ) 
represent two half-order differentials, then higi dz extends meromorphically to 
fl as a meromorphic differential (via extension by h 2 g 2 dz on the reflected side). 
We let hgdz denote the extension of this form to the double. Hence, if 7 is a 
curve on fl and h = (hi, h 2 ) and g = (gi,g 2 ) are half-order differentials, we may 
define the pairing 



B(h,g) = 




as long as the extension hg dz to fl does not have poles along the curve. 

Notice that identity ( 15.11) is the same as the critical identity ( 13. ip for the Szego 
and Garabedian kernels. Results in [7] (Lemma 6.1) show that all meromorphic 
half-order differentials without poles on the boundary arise as linear combinations 
of the Szego and Garabedian kernels and their derivatives in the second variable. 
When there are no poles inside the domain, a meromorphic half-order differential 
must be a linear combination of just the Szego kernel and its derivatives in the 
second variable on the domain side. This is the key connection between the Szego 
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kernel and the results from [T5] that allow us to use Szego coordinates to prove 
the density of double quadrature domains here. 

The idea of the proof is to construct a half-order differential h = (hi, h 2 ) with 
the property that hi is in the Szego span and is close to one in A°°(Q), and h 2 dz 
is a meromorphic one-form on Q with a single pole at a point in the reflected 
side of Q in Q with the property that h 2 dz = for j = 1, . . . , 2p. Then 

h 2 dz = df where / is a meromorphic function on f2. If hi is sufficiently close to 
one, then the condition /' = h 2 implies that / is a biholomorphic mapping of f2 
onto a domain which is C°° close to Q. Now Theorem 11.11 yields that f(£l) is a 
boundary arc length quadrature domain. Finally, by [Hj, the fact that / extends 
meromorphically to the double yields that /(fi) is an area quadrature domain. 
In this way, we will obtain a double quadrature domain that is C°° close to Q. 

We shall need the following special Mergelyan approximation lemma. Suppose 
a G U. Recall that S ao denotes the complex linear span of functions of z from 
{S m (z,a ) :m = 0,1,2,...}. 

Lemma 5.1. Given a holomorphic function G in A°°(Q) and continuous func- 
tions (fj on Tj that satisfy the compatibility condition 

G(w)T(w) = <fj(w) 

at the endpoints w of each curve Tj (which fall on bfl), it is possible to find a 
function 

N 

a i z ) = c mS m (z,a ) 

m=0 

in S ao such that a is as close as desired to G in A°°(Q) and so that 

N 

A(z) = -i^2~c^L m (z,a ) 

m=0 

is as close to tpj on each Tj in the uniform topology as desired. (Note that identity 
(COP shows that (a, A) is a half- order differential.) 

We will now prove the density of double quadrature domains, assuming the 
truth of the approximation lemma. Afterwards, we will prove Lemma 15.11 

Let K denote the compact set Vt U (u^ =1 7 J+p ) in Vt. We may consider non- 
meromorphic half-order differentials defined on K where g = (g\,g2) is such that 
gi is in C°°(Q) and the g% is a continuous function on the union of the curves Tj 
and satisfies the compatibility condition 

g x {w)T(w) = g 2 (w) 

at the endpoints w of each curve Tj. For two such half-order differentials g and 
h, we may define the pairing 



B k (g,h) = 
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for k — 1, . . . , 2p as described above in the meromorphic case. 

We now define 2p+l such half-order differentials g m = (g™, g™), m — 0, . . . , 2p, 
on K satisfying 

(i) 5 ? = lon O, 

(ii) B fc (c/ , #°) = for all k = 1, ... , 2p, and 

(iii) B k (g°, g m ) = 5 km for k, m = 1, . . . , 2p, 

We start by setting (7° = 1 on Q (where g® is the first function associated to the 
half-order differential g°) as condition fli]) dictates. Notice that (JTTJ) is therefore 
automatically satisfied for k = 1, . . . ,p. It is now an easy matter to define g® on 
the curves Tj so that the compatibility conditions hold and condition (jn]) holds 
for k — p+1, ■ ■ ■ , 2p as well. Next, set g™ = R m on Q, where R m are the rational 
functions of the previous section for m — 1, . . .,p. Note that condition flm|) now 
holds for k, m — 1, . . . ,p. It is another simple matter to define g™ on the curves 
Tj so that ( lull) also holds for k = p + 1, . . . , 2p and m = 1, . . . ,p. Finally, g m are 
chosen for m = p+1, . . . , 2p so that = on and g™ are defined on the curves 
Tj so that flm|) holds for k —p+1, ... , 2p. Note that flmj) holds automatically 
for = 1, . . . ,p. 

Lemma T5. II implies that there are half-order differentials h m = (/i™, h™) which 
are close to the g m in the sense that h™ is close to g™ in A°°(Q) for each m and 
h™ is uniformly close to g™ on the curves Tj. The half-order differential that we 
seek will be given by 

2p 

h = h° + J2 t m h m , 

m=l 

where the t m are complex numbers close to zero. (Note that we may take linear 
combinations of half-order differentials in the obvious way by summing the corre- 
sponding functions in the pairs.) We will choose the t m to make B k {h, h) = for 
k = 1, . . . , 2p so that h 2 dz is period free, and is therefore an exact meromorphic 
differential df. (Note that the lone pole of h 2 dz would necessarily have residue 
equal to zero.) Let g = g° + Y^m=x tmg m ■ Define complex coefficients via 

a kij = Bkih* + g\h j - g j ) 
b kj = 2{B k {gi,ti ) -g°) + B k {ti 3 ,y - g*)) 
c k = B k (h° + g ,h°-g ) 
A kij = B^g^gi) 

Note that A ki j is fixed, while the others depend on the h? and can be made 
arbitrarily small by improving the approximation sufficiently. Expand B k (h, h) = 
B k {h — g + g,h — g + g) to obtain 

B k (h, h) = 2t k + c k + ^2 AkijUtj + ^2 

i,j i,j j 

We obtain a system by setting these 2p periods equal to zero. Note that this 
system has the unique zero solution if the coefficients that depend on the h? 
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are taken to be zero. It is now a standard application of the implicit function 
theorem to deduce that it is possible to choose tj close to zero to make all these 
periods vanish. This completes the proof, assuming Lemma 15.11 We now turn 
to proving the lemma. 

Note that Theorem II .21 allows us to simplify the assumptions in Lemma EH] by 
subtracting off an element of the Szego span to be able to assume that G is zero 
on Q. Given the point ao in the simply connected domain U — Q — U p =l Tj, our 
task then is to find a function X(z) which is a linear combination of functions of 
the form L m (z,ao) which is uniformly close to the functions tpj on each Tj and 
which is close in C°°(bQ) to zero. This turns out to be remarkably easy. The 
first step in the argument is to use the classical Mergelyan Theorem (see Rudin 
[IS] . Chapter 20, Exercise 1) and Runge's Theorem to approximate our functions 
as follows. For e > 0, let V e denote the points in C that are less than or equal 
to a distance e from bQ, i.e., the closure of a small collared neighborhood of the 
boundary in C. Let K e = V e U (u P j =1 Tj). We first use Mergelyan's Theorem 
to obtain a rational function R(z) with poles in C — K e that is unformly close 
to zero in V e , and because the functions (pj vanish on bQ, and because we may 
shrink e, we may also arrange for R(z) to be uniformly close to the (pj on each 
Tj. Cauchy's Estimates reveal that uniform convergence on the interior of V e 
implies C°° convergence on bQ. Hence, we may obtain a rational function that 
is close to zero in C°°(bQ) and uniformly close to (pj on each Tj. Next, we may 
slide the poles of R(z) that fall in U to the point ao as is done in many standard 
proofs of Runge's Theorem (see Stein [221 P- 63]) to obtain a rational function 
R(z) with poles at a and finitely many other points outside of Q which is as 
close in C°°(bQ) to zero as desired and as close in the uniform topology to (fj on 
each Tj as desired. 

The Garabedian kernel L(z,w) is such that 

L(z, w) = r + £(z, w), 

2n(z — w) 

where £(z, w) is holomorphic in z and w and is in C°°(Q x Q) (see [5j p. 102]). We 
will therefore be able to use L(z, w) much like a Cauchy kernel in what follows. 
The Residue Theorem yields that, for z G Q — {a }, 

1 f 

- / L(z,w)R(w)dw 
i Jwebn 

is equal to R(z) plus 27r times the residue of L(z,w)R(w) at a , which is a fixed 
finite linear combination of functions of the form L m (z, ao). If we show that the 
integral tends to zero in z in A°°(Q), then we will have completed the proof. 
This turns out to be an easy consequence of the fact that the Szego projection 
is a continuous linear operator in C°°(Q) (see p. 13]). Indeed, identity (13. II) 
and the fact that L(z,w) = —L(w,z) yields that the integral is equal to 




S(z, w)R(w)ds, 
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which is the Szego projection evaluated at z of the function which is the restric- 
tion of R(w) to the boundary. Since R is close to zero in C°°(bQ), it follows that 
the function given by the integral is close to zero in A°°(bQ). The proof of the 
lemma is finished. 

6. Properties of double quadrature domains 

Suppose Q is a double quadrature domain in the plane. Because Q is an area 
quadrature domain, it has a boundary given by an algebraic curve, and because 
it is a boundary arc length quadrature domain, that curve consists of C°° smooth 
real analytic curves (see [TJ [HI [15]). Because Q is an area quadrature domain, the 
complex polynomials belong to the Bergman span (see [H]). Also, the Schwarz 
function S(z) exists and is meromorphic on Q, extends analytically past the 
boundary, and satisfies S(z) = z on the boundary (see [H [33]). This identity 
shows that the two functions z and S(z) extend meromorphically to the double of 
Q. It was proved in [H] that the extensions form a primitive pair for the field of 
meromorphic functions on the double, i.e., all the meromorphic functions on the 
double are given as rational combinations of the two. It is shown in [8] (see also 
[9]) that the Bergman kernel K(z,w) associated to an area quadrature domain 
is a rational function of z, S(z), and w, and S(w). Consequently, K(z,w) is a 
rational function of z, z, w, and w when z and w are restricted to the boundary. 
It was also shown in [8j [9] that the Bergman kernel is a rational combination of 
two Ahlfors maps associated to two generic points in the domain. Aharonov and 
Shapiro pQ showed that the Ahlfors maps are algebraic. 

Because Q is an arc length quadrature domain, there is a meromorphic function 
H(z) on Q which extends analytically past the boundary and such that H(z) is 
equal to the complex unit tangent vector function T(z) on the boundary (see 
[T5l [2]). The identity relating the Szego and Garabedian kernels can be used 
with this fact to see that both the Szego and the Garabedian kernels extend to 
the double in the first variable when the second variable is held fixed. Indeed, 
the identity 

S(z,a) = -L(z,a)T(z) for z G bQ 
i 

shows that L(z, a) = iS(z, a)/ H(z) on the boundary, and this shows that L(z, a) 
extends meromorphically to the double in z. Similarly, S(z, a) extends. These 
identities also reveal that T{z) is equal to a function G(z) on the boundary 
where G(z) is meromorphic on Q and extends analytically past the boundary. 
Hence, T(z) is the restriction to the boundary of a meromorphic function on the 
double of Q that has no singularities on the boundary of Q. These arguments 
are reversible and we may state the following theorem. 

Theorem 6.1. A bounded domain with a piecewise C 1 smooth boundary is a 
boundary arc length quadrature domain if and only if the Szego kernel S(z, w) 
associated to the domain extends meromorphically to the double as a function of 
z for each fixed w in the domain. This is the case if and only if the complex 
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unit tangent vector function is the restriction to the boundary of a meromorphic 
function that has no singularities on the boundary offl. 

The same theorem holds with the Garabedian kernel in place of the Szego 
kernel. (The part of the theorem about T was proved in [15]. The equivalence 
of the extendibility of the Szego kernel is a rather direct consequence of the 
extendibility of T and might very well have been noted earlier.) 

If is a double quadrature domain, then we have shown that T(z) is the 
restriction to the boundary of a meromorphic function on the double, and we 
also know that such functions are generated by z and S(z). Hence, it follows that 
T(z) is the restriction of a rational combination of z and S(z) to the boundary, 
i.e., T(z) is a rational function of z and z. 

When Theorem 16.11 is combined with results from [B] , it follows that the Szego 
kernel S(z, w) associated to a double quadrature domain is a rational combination 
of z, S(z), and w, and S(w). Consequently, S(z, w) is a rational function of z, 
z, w, and w when z and w are restricted to the boundary. 

The Kerzman- Stein kernel is given by 



for z and w in the boundary of f2, and it follows that the Kerzman-Stein kernel 
associated to a double quadrature domain is a rational function of z, z, w, and 
w. 

We shall now show that the complex polynomials are in the Szego span associ- 
ated to a double quadrature domain. Let H (z) denote the meromorphic function 
that is equal to T(z) on the boundary. Suppose h G A°°(Q). Since z = S(z) on 
the boundary, we may write 



and this last integral is equal to a finite linear combination of values of h and 
its derivatives at finitely many points in Q by the Residue Theorem. Hence, z n 
has the same effect when paired with h that a certain element in the Szego span 
would have. Since A°°(Q) is dense in the Hardy space, it follows that z n is equal 
to the element in the Szego span. 

The converse of this last result is true, namely, that if the complex polynomials 
are in the Szego span, then the domain is a double quadrature domain. To see 
this, note that z = z/1. If both z and 1 are in the Szego span, then z is a 
quotient of elements in the Szego span, and as such, it extends to the double 
as a meromorphic function (since the relationship between the Szego and the 
Garabedian kernels reveals that quotients of functions in the Szego span would 
be equal to quotients of conjugates of functions in the "Garabedian span" on the 
boundary). The condition that z extends to the double is the classic condition 
that is equivalent to the domain being an area quadrature domain. It is easy to 
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see that a smooth domain is a boundary arc length quadrature domain if and 
only if the constant function 1 is in the Szego span. Hence, we have proved the 
following theorem. 

Theorem 6.2. A bounded domain with a piecewise C 1 smooth boundary is a 
double quadrature domain if and only if the complex polynomials belong to the 
Szego span associated to the domain. 

Theorem 16.21 is an appealing analogue of Theorem 1.2 from [TT] which states 
that a bounded domain is an area quadrature domain if and only if the complex 
polynomials belong to the Bergman span. 

The Poisson kernel of a smooth simply connected domain is given by 



where z e Q, w G bQ, and a is a fixed point in Q (see P, p. 1367]). Hence, when 
Q is a simply connected double quadrature domain, p(z, w) is the real part of a 
function that is rational in z and S(z), and w and w (just like in the unit disc, 
where S(z) = 1/z). 

The rest of the properties mentioned in §1 involving the Riemann map / of 
a simply connected domain follow from the remarks above and the fact that 
/ extends meromorphically to the double and the extension generates the field 
of meromorphic functions on the double, i.e., every meromorphic function on 
the double is a rational function of the extension of the Riemann map. Hence, 
for example, on a quadrature domain with respect to area, the functions z and 
S(z) extend to the double, and so they are rational functions of f(z). Since, 
f{w) = 1/ f(w) on the boundary, the restriction to the boundary of a rational 
function of f(w) and f(w) is equal to a rational function of f(w) alone. Hence, 
it follows from the remarks above about the Poisson kernel is the real part of a 
rational function of f(z) and f(w). 

7. The Schwarz function and proper holomorphic mappings to the 



We conclude this paper by noting a relationship between proper holomorphic 
mappings of a domain to the unit disc and the existence of the Schwarz function, 
which is the hallmark of a bounded quadrature domain with respect to area. We 
go on to reveal a further relationship between proper holomorphic mappings and 
the extension of the unit tangent function to the double, which characterizes a 
bounded quadrature domain with respect to area. 

Assume that Q is a (bounded) area quadrature domain. For the moment, 
assume that zero does not belong to Q. As remarked earlier, the Schwarz function 
S(z) exists and is meromorphic in Q. Let {aj}^ =1 denote the poles of S(z) in 
fl, and suppose aj is a pole of S(z) of order rij. Let fj denote the Ahlfors 
map associated to aj, which is a proper holomorphic mapping of Q onto the 
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unit disc such that aj is a simple zero of fj. (The Ahlfors map associated to a 
point a in Q is the unique solution to the extremal problem to maximize h'(a) 
under the conditions that h maps Q into the unit disc and h'(a) is real. It maps a 
bounded n-connected domain onto the unit disc as an n-to-one branched covering 
map.) We note here that, since the boundary of Q is piecewise smooth, proper 
holomorphic mappings to the unit disc extend continuously to the boundary. 
Conversely, a holomorphic mapping to the unit disc that is continuous up to the 
boundary is proper if and only if it is non-constant and it maps the boundary to 
the unit circle. Notice that, since = 1 and z = S(z) on the boundary, the 

function 



has removable singularities at each a J; is continuous up to the boundary, and has 
unit modulus on the boundary. If Q is multiply connected, then $ must have a 
zero in Q since the Schwarz function must have a zero in the domain (see [2]). 
Hence $ is a proper holomorphic mapping of Q to the unit disc. It follows that 
S(z) is z times the quotient of two proper holomorphic mappings of Q to the 
unit disc. If the point zero were in the domain, we could add the Ahlfors map 
fo to the product and obtain the same result. If Q is simply connected, then we 
can define 



where fj is the Riemann map associated to a point aj where S(z) has a pole of 
order rij, and is the Riemann map associated to a point bk where S(z) has 
a zero of order m^. In this case, \I/ has constant modulus one on the boundary 
and no zeros in Q. Hence it is a unimodular constant and we deduce that S(z) is 
equal to a unimodular constant times z times a quotient of products of Riemann 
maps. Since proper maps are all given by unimodular constants times products 
of Riemann maps in this setting, the same result follows, i.e., that S(z) is z times 
the quotient of two proper holomorphic mappings of Q to the unit disc. 

If the domain Q is bounded by finitely many non-intersecting Jordan curves, 
a simple converse follows. Indeed, since a proper holomorphic mapping / to the 
unit disc would be continuous up to the boundary and have unit modulus on the 
boundary in this case, the identity f(z) = 1/ f(z), which holds on the boundary, 
reveals that / extends to the double. If z is equal to z times the quotient of two 
proper holomorphic mappings to the unit disc on the boundary, then z is seen 
to extend to the double, and it follows that the domain is an area quadrature 
domain. The theorem with its converse reads as follows. 

Theorem 7.1. A bounded domain in the plane bounded by finitely many non- 
intersecting Jordan curves is an area quadrature domain if and only if, on the 
boundary, the function z/z is equal to the boundary values of a quotient of proper 
holomorphic mappings of the domain onto the unit disc. 
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Note that if z/z — /1//2 on the boundary, then the Schwarz function is given 
by S(z) = zf 1 (z)/f 2 (z). 

If a domain has piecewise C 1 smooth boundary, it is a boundary arc length 
quadrature domain if and only if the unimodular function T(z) extends to the 
double as a meromorphic function. Similar reasoning to the argument above 
with T(z) in place of S(z)/z yields the following result. 

Theorem 7.2. A finitely connected bounded domain with piecewise C 1 smooth 
boundary is a boundary arc length quadrature domain if and only if the complex 
unit tangent vector function T(z) is equal to a quotient of proper holomorphic 
mappings of the domain onto the unit disc. 

Since an area quadrature domain has piecewise C l smooth boundary, we may 
also state this last theorem. 

Theorem 7.3. A bounded domain in the plane bounded by finitely many non- 
intersecting Jordan curves is a double quadrature domain if and only if, on the 
boundary, the function z/z is equal to the boundary values of a quotient of proper 
holomorphic mappings of the domain onto the unit disc, and the complex unit 
tangent vector function T(z) is also equal to the boundary values of a quotient of 
proper holomorphic mappings of the domain onto the unit disc. 

Finally, we remark that the semi-group of all proper holomorphic mappings 
of a finitely connected domain to the unit disc has been described completely in 
US- 
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